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ABSGEACT 



The standard deviation computer is a small electronic analog device -which 
is designed to calculate the standard deviation of a voltage signal with less 
than 5^ error. The output is a continuous voltage signal proportional to the 
standard deviation which is generated with a computation time equal to the 
sampling period. Mathematically^ this instrument computes the mean absolute 
deviation which is related to the standa^i'd deviation "by a constant factor • It 
accepts signals in the frequency range of from d-c to one kilocycle over a 
wide range of amplitudes and computes with sampling times up to 560O seconds 
or more. In field tests, the computer operated successfully on the output of 
a wind vane . 
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INTRODUCTION 



Progress in the study of atmospheric turbulence has been hindered by the 
lack of satisfactory means of measuring directly the significant turbulent 
characteristics of the wind g.t any given time and place. Earlier methods^ such 
as those used by Scrase and Best^^ required analysis of the traces made by a 
sensitive anemometer^ wind vane, or bivane. In one method the "mean fluctua- 
tion ratio," which is the mean of the absolute values of the eddy velocity ex- 
pressed as a percentage of the mean wind speed for a given period, is deter- 
mined by a time-consuming analysis of wind-speed chart traces. The "gustiness" 
or "intensity of turbulence," which is the ratio of the root -me an -square value 
of the oscillations of the wind speed to the mean wind speed, is also cumber- 
some to evaluate . The gustiness or intensity may also be determined approx- 

p 
imately from fluctuating wind direction traces, as shown by Best, but this is 

also time-consuming and difficult to establish on a routine basis so that the 
evaluation of records may be done by nonscientific personnel. 

The first instrinnent designed specifically to give a direct measurement 
of atmospheric turbulence was the bridled-cup turbulence integrator or bridled - 
cup gust accelerometer developed by Hewson and Gill.^^^ This relatively simple 
instrument gives the total acceleration, without regard to sign, occurring in 
the fluctuating wind speed during any specified time interval. For example, 
the instrument may show a total wind acceleration of 800 m sec^ during a pe- 
riod of 50 min, and this value may be obtained directly from the chart by un- 
skilled workers without any exercise of scientific judgment on their part. The 
instrument thus gives a measure of a fundamental property of turbulent flow, 
its acceleration, and of the turbulent diffusion associated with it. Its gen- 
eral usefulness has been limited, however, by the fact that the theories of 
atmospheric turbulence and diffusion have not yet developed substantially in 
terms of the accelerations of turbulent flow, although a beginning in this dir- 
ection has been made by Cramer .5 

There has been a recent trend toward the specification of the turbulent 
characteristics of the atmosphere in terms of the statistical properties, such 
as the standard deviation, of the wind speed and direction. !]liis approach was 
developed by Hewson, Gill, Cramer, and Record,^ and used by Friedman, ' Falk 
et al.,^ and Record and Cramer-^ in field studies of atmospheric diffusion. 
Further development and comparison with field results have been undertaken by 
Hay and Pasquill; and Panofsky and McCormick-*--^ have suggested a semi -empirical 
expression for the variation with height of the standard deviation of the ver- 
tical component of the wind velocity. 

A direct method of computing an estimate of the standard deviations of the 
wind direction from some arbitrarily fixed bearing has been presented by Jones 
and Pasquill.-^^ The experimental system reported by them was designed to com- 



pute the statistical mean of the standard deviation of wind-direction fluctua- 
tions -when averaged over 5 seconds and sampled over 5 minutes* 

Study of the fundamental properties of atmospheric turbulence and diffu- 
sion vould "be facilitated by a similar device capable of calculating the stand- 
ard deviation -with less than 5^ error and reliable enough for continuous opera- 
tion. This accuracy requirement is imposed to match the accuracy of the sensor 
and so that the data need not be recomputed* The computer should provide sev- 
eral pre-set values of the sampling period up to 560O seconds. Since the com- 
puter is to perform mathematical operations on the signal, it should be pos- 
sible to convert it readily to other uses, such as computing the correlation 
function of two signals . On-site operation demands that the computer be com- 
pact and rugged with simple but effective controls . 

Tbese requirements suggest the use of an electronic analog circuit which 
would consist of high-quality passive components, e.g., resistors and capacitors, 
and of low-cost, commercially available operational amplifiers. Using these 
components, precise mathematical operations can be performed and this type of 
computer will accept any direct analog voltage signal and will drive any con- 
ventional recorder or display device. 



2. MATHEMATICAL DESIGN CRITERIA 



The mean and standard deviation of a random signal, x(t), over a fixed 
period of time are given "by 
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= \J- x(t) dt 



(1) 



s = < 
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where 



x(t) = input signal 

t = time 

T = sampling period 

X = mean over period t 

s = standard devieition over period t. 



Another measure of the dispersion which Thrill be useful is the mean absolute 
deviation given by 



E 



P+1 



(3) 



where E is the mean absolute deviation over the period t. If x(t) is normally 
distributed, the standard deviation and mean absolute are related by 



s = 1.25 E 



i}^) 



This relation may also hold for some non-normal distributions but if {h) does 
not hold, the distribution is not normal. ^ 



The above foarmulae may be modified for continuous computation by changing 
the limits of integration to t - (t/2) and t -h (t/2) instead of -t/2 and +t/2, 
respectively. At this point, a difficulty arises because no physical system can 



perform Integration in future time as required by the limit t + (t/2). This 
difficulty can he overcome hy delaying the computation by t/2 seconds so that 
the mean is given by 



x(t - t/2) = i 



x(t) dt 



(5) 
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The expressions for the standard deviation and the mean absolute deviation are 
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E(t - t) = 



t, 



|x(t - t/2) - x(t - t/2) I dt 



t-T 



(7) 



Equations (6) and (7) require comparison of the input with the mean of the in-| 
put computed over a period symmetrical ahout the input so that if the mean is 
given by (5)^ the instantaneous value to he compared with it is x(t - t/2) 
rather than x(t). Delay of the input by t/2 seconds for large t requires ex- 
pensive equipment not compatible with the demand for low initial cost. What 
is needed is an approximation to x(t - t/2) - x(t - t/2) which is sufficiently 
accurate but which can be implemented with a simple electrical analog. If the 
approximation can be stated as a ratio of two polynomials in the Laplace trans- 
form domain, in the form 
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where S is the transform variable, then there exists a direct electrical analog 
which can be built up with relatively inexpensive components. 



The Laplace transform of x(t - t/2) - x(t - t/2) is 



x(t - t/2) - x(t - t/2) 
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Q3ie Infinite series expression for (9) is 



n+1 n+2 



X(S - t/2) - X(S - T/2') ^ ^ (-1) (tS) 
X(S) n=0 
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which can be approximated by the ratio of two polynomials in S of the form of 
(8). For a second-order polynomial^ the approximation is 



X(S ■ t/2) - X(S - t/2) ^ 0.067t^ S^ 

X(S) 0.067t^ S^ + O.5OTS + 1 



(11) 



Oliis was obtained by graphical fitting. The result may be seen in Fig. 1. 

The exact solution^ needed for comparison^ is more readily obtained by 
direct integration than by inversion of the transform given in (9). For a 
sinusoidal input of unity amplitude^ 



x(t) 






sin a)t dt = 
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and 



X 



(t - T/2) - x(t - t/2) = 



(JDT - 2 sin 2 ^^ 



CDT 



sin a)(t - t/2) (15) 



The term in brackets in (13) is the amplitude of the variations about the mean 
and the phase shift is just t/2 seconds. 

Equation (11 ) may also be put in the amplitude form similar to (13) by 
converting to a Fourier transform. This is done by replacing S by jo). Then 
the amplitude term is given by the real part: 
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(Hie correct amplitude function for x(t - t/2) - 5?(t - t/2) given by the term 
in brackets in Eq.. (I5) is plotted in Fig. 1 as a, function of odt along with 
the corresponding amplitude function for the approximation given by Eq.. (lla) . 

Ihe second integration to be performed, the one indicated in Eq.s . (6) and 
(7)^ must also be put in the form of (8). If this operation is 



Yi dt 
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then the Laplace transform is 
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An approximation to (13) can be developed from the series expansion: 
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n=0 (n + 1): 
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The second-order approximation to (15) is 



12 



T^ S^ + 6 tS + 12 



(16) 



The desired amplitude function for a unity sinusoidal input is given by the 
term in brackets in Eq. (12). This is plotted in Fig. 2 as a function of cm 
along with the real part of the Fourier transform of (l6), which is 
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The mean absolute deviation for the vmity amplitude sinusoidal input can 
be found by combining Eqs( ,(7) an4- (13)^ 
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Fig. 2. Plot of the Integration of sin cot over sampling period 
T with the second-order approximation. 



8 



E 



(t - t) = COT - 2 sin 1/2 (OT /" 1^.^ ^(^ _ ^/2^| 

CDT^ ^t-T 



dt 



(17) 



The absolute value term can be expressed as a Fourier series 
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This expression will be used to compare vith the frequency-response data of 
the computer . 



COMPUTING CIRCUITS 



The computing circuits required for generation of the mean absolute devia- 
tion and the standard deviation were set up on a small ^ general -purpose analog 
computer to demonstrate the feasibility of the system* This is a practical test 
since this computer has the same quality components and amplifiers as the pro- 
posed special purpose computer. The circuit diagrams are shown in Figs. 5 
and k. All the indicated resistors and csipacitors are 1^ quality and the multi- 
pliers are the diode quarter -square type. The circuit diagram of Fig. 5 ttas 
been divided into two major mathematical operations by the dashed lines. The 
operations are the generation of |x(t - t/2) - x(t - t/2) | in the top half of 
Fig. 5^ 9J^^ "t^^^ subsequent integration required for computation of E at the 
bottom^ The coefficients are determined by three circuit components^ a volt- 
age divider a, p^ d, or y followed by an E.C combination as follows: 



3>86 Ri Ci ^ IM Rg Ci ^^^^ 

a p 



from Eq. (12)^ and 



3-^7 Rs Cs 6 R4 Cs ,. 

T = [ = (20) 



from Eq. (16) . Equations (19) and (20) must be satisfied for the concept of 
sampling time to have meaning. 

Given that the computer can handle a voltage range from -100 to +100 volts, 
then, for optimum use of the computer, the voltage input x(t) should have ex- 
treme values of ±100 volts. If this is not practical, the input resistor R^ 
shown in Figs. 3 and h provides an adjustable gain. If the voltage extremes of 
x(t) are, say, ±10 volts,- then Rg should be 0.1 megohm which makes the effec- 
tive voltage range ±100 volts since the input gain is l/Rg. 

The relation between the voltage at any point in the circuit and the mag- 
nitude of the corresponding physical variable is determined by the amplitude 
scaling factor A which is defined by the relation 



X(t) = A x(t) (21) 



where x(t) is the physical quantity and X(t) is the corresponding voltage. This 
relation holds at every point in the circuit in Figs . 5 and k . Ordinarily, the 

11 



— ^AA»- 



Ra 



> 



/ 

-(X-X) 






r-'VsAr 



^'{^^^^^^'hz® 



+(X-X) 



-(a)-WSr 



Rz 

■AAA* — ' 



He 



Ci 



^ 



i-WAr 



J'-^^l 25 Cz Cz 

-AA^ pHe^ r-H(— 



IX-XI 



1.25 E 




HIGH-GAIN d-c AMPLIFIER 



VOLTAGE DIVIDER 



Fig. 5- Computing circuit for generation of 1.25 E vhich is an 
estimate of s. All resistor and capacitor values are in units 
of megohms and microfarads^ respectively. 
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squaring operation changes this relation but the circuit points affected have 
been relabeled so that (21) holds throughout. 

Ihe circuit of Fig. 4, for computing the standard deviation, is similar 
to that of Fig. 3 but has a squaring device instead of an absolute value net- 
work and another multiplier set up to perform the square root operation at the 
end of the circuit . 

To compute the standard deviation, 8 amplifiers and 2 multipliers are 
needed but only 7 amplifiers are required for the mean absolute deviation. 
Simply on the basis of the amount of equipment required, it is least expensive 
to compute the mean absolute deviation. 

It is possible to design a circuit to perform the same functions "with 
fewer amplifiers. A ^-amplifier circuit for computing the mean absolute devia- 
tion is shown in Fig. 5 "v^here the transfer function of the first amplifier is 



2 

fo _ Ci Cg Ri Ra S , . 

^i Ri R2 Ca C3 S^ + Ri(Ci + C2 + 03)8 + 1 



and Ci = C3 for unity gain. To determine parameter values, Eq. (22) must be 
matched with Eq. (ll). The transfer function of the last amplifier is 



fo ^ R5 1 . . 

^i R3 R4 R5 C4 C5 S2 + ^ (R4 R5 + R3 R4 + R3 R5)S + 1 

and R5 = 1.25 R3 for a gain of 1.25 . Equation (25) must be matched with Eq.(l6) 
This circuit has the disadvantage that it does not work well for values of the 
sampling period greater than about 10 seconds. An amplifier with an output 
capacity of 20 milliamp at ±100 volts could not drive the complex feedback cir- 
cuits for T greater than 12 seconds. Also, the magnitude of the resistors and 
capacitors required is much greater than for the 7-QJ3iplifie^ circuit. For com- 
parison, the 7-amplifie3? circuit, using the same amplifiers, can operate with 
a sampling period of 36OO seconds and with reasonable component values. For 
example, a t of 360O seconds can be obtained with 10-megohm resistors, 10- 
microfarad capacitors, and voltage divider settings between O.O9 and 0.17' 
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h. EVALUAHON OF THE COMPUTER 



Evaluation of the computer was directed toward justifying the use of the 
mean alDSolute deviation as an approximation to the standard deviation. The 
computer performance was evaluated with respect to Eq. (T) an<3. an attempt was 
made to evaluate Eq. (4). The procedure involved both laboratory and field 
tests .. 



A. LABORATORY TESTS 

Since any input can "be made up of a combination of sine waves of various 
frequencies, it is valid to study the frequency response of the computer. A 
sine wave generator was used as the input;, and to facilitate the work, a sam- 
pling period of 2 seconds was used. The result was plotted in Fig. 6 in terms 
of CDT so that the data apply to any sampling period. The output of the com- 
puter was compared with the theoretical response given by Eq. (l8), which is 
the integrated form of Eq. (7)^ where the input is a unity magnitude sine wave. 

The mean absolute deviation of a sine wave consists of a steady component 
plus a component which oscillates in time. This is indicated in Fig. 6 by 
plotting the mean or steady value and the bounds of the oscillating components . 

For 0.1 ^ CDT ^ 100, the computer error relative to Eqs . (ll) and (l6) de- 
pended upon the amplitude of the output so that the error was less than 5^ 
at the 5 -volt level, less than 3?^ s-t the 10-volt level, and less than 2^ at 
the 50-volt level. Error depending on the amplitude level is a phenomenon com- 
mon to all such electrical computing devices • 



B. FIELD TEST 

As a further demonstration of the computer perfoarmance, it was given a 
field test. A potentiometer -type wind van.e was mounted on the roof of a build- 
ing. Voltages applied to the potentiometer were ±10 volts from the computer, 
and the potentiometer output was fed directly to the computer with a gain of 
10 so that ±l8o degrees of vane rotation equals ±100 volts • Using the circuits 
shown in Figs* 3 and 4, the standard deviation and the mean absolute deviation 
were computed. The sampling times used were 2k and 6o seconds. The amplitude 
scale factor A, was I/I.8. 

This was a fairly severe test as the wind vane was in a region of mechan- 
ical turbulence, although complete revolutions of the vane were uncommon. The 
standard deviation was in the range of from I5 to 20 degrees. No provision 
was made to keep the vane from rotating the potentiometer arm across the ends 
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of the potentiometer winding "which causes the amplified computer input to jump 
from +100 volts to -100 volts or vice versa • To the computer^ this looks like 
part of the signal upon which it is to operate^ and this causes a large per- 
turbation in "both measures of the deviation for about 5't seconds . 

Some of the results were recorded with a chart speed of 10 millimeters 
per second, which is fast enough to peratit good resolution of the input signal. 
From these, E was hand -calculated to check the computer. 

This calculation, over a 45 -second portion of the record where t = 2^ 
seconds, showed an average k^ error in the mean absolute deviation. Instantan- 
eous values of the error were as great as 8^, but this is at least partly 
attributable to abstracting error. 

It was stated above, Eq. (4), that E and s should be related by the form 
factor 1.25 when the input is normally distributed. Having computed both E 
and s, this can be readily checked. To do this, E and s/E were determined at 
a number of points in an interval of the record. The fonm factor and the mean 
absolute deviation were averaged for the points to eliminate random measure- 
ment errors as shown in Table I . 



TABLE I 



FORM FACTOR FOR SEVERAL VALUES OF TEE SAMPLING TIME 



Sampling 


Number 


Length of 


" 


time, 


of 


record, 


E, 


sec 


points 


sec 


deg 


2k 


20 


20 


li^.i^ 


24 


20 


20 


Ik.^ 


60 


20 


ko 


12.7 


60 


20 


ko 


li^.6 



Form 
Factor 



1.21 
1.21 
1.23 
1.24 
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SUMMARY AND CONCLUSIONS 



Hiere were two key steps involved in this method of computing the standard 
deviation. The first was to find a mathematical model which is an adequate 
approximation to the defining model and which could he readily implemented as 
an electrical analog. The second was to find this electrical analog which has 
all the desired properties^ ^-g*^ accuracy^ simplicity^ and reliability. 

The defining model was expressed in Eqso (4) and (7) and the approximate 
model was partially stated in Eq.s« (11 ) and (16). The complete statement is 
that the input is operated on "by (ll); the absolute value is taken; this is 
operated on "by (l6)j and the result is mu]-tiplied "by 1.25 to obtain the stand- 
ard deviation. 

Computer circuits have been shown (FigSo 5 and k) for obtaining the stand- 
ard deviation by two methods . The first used the concept of the mean absolute 
deviation and its relation to the standard deviation^ while the second circuit 
obtained the standard deviation directly with the variance as a byproduct. The 
second circuit has the advantage of being somewhat more exacts but it is also 
more expensive to implement. The decision to use the first circuit was made 
solely on the basis of initial cost. As another measure of economy^ the cir- 
cuit of Fig. 5 can be used^ provided that only small values of the sampling 
time are required. 

This computer can accept any direct analog voltage signal at any ampli- 
tude level. It can attenuate or amplify that signal to its own optimi;ijn range. 
The upper limit in frequency of input is the frequency at which the amplifier 
band width interferes^ so that the upper limit should be at about one kilo- 
cycle. The maximum sampling time that can be used is ultimately set by the 
capacity of the amplifiers to drive the computing networks and by the diffi- 
culty in avoiding leakage paths around very large resistors . With amplifiers 
having a forward gain of 200,000 and an output capacity of 20 milliamp at ±100 
volts, the accuracy will be unimpaired with sampling periods up to 560O sec- 
onds but will begin to fall off for larger sampling periods . 

The primary limit on accuracy is the quality of the mathematical approxi- 
mations used. The computing circuits were implemented with 1^ quality re- 
sistors and capacitors, but if more accuracy in the computer were needed, one 
could easily use 0.05^ components. 

Stability of operation is an inherent feature, even though the d-c ampli- 
fiers used are not chopper-stabilized • The computer can operate unattended, 
with its accuracy unimpaired by amplifier drift, for periods of up to one week. 
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The applications for this computer are n-umerous . In addition to computing 
the standard deviation for air -pollution studies, it could "be used as an aid 
in the controlled emission of plant effluents, the study of time averaging of 
the standard deviation of turlDulent eddies versus spatial averages, quantifying 
the turbulent structure of air flow around louildings, and the study of the de- 
cay of turbulence with height* 
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